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ABSTRACT 

We show that collision-induced absorption allows molecular hydrogen to act as an inconden- 
sible greenhouse gas, and that bars or tens of bars of primordial ti2-He mixtures can maintain 
surface temperatures above the freezing point of water well beyond the "classical" habitable zone 
defined for CO2 greenhouse atmospheres. Using a 1-D radiative-convective model we find that 
40 bars of pure H2 on a 3 Earth-mass planet can maintain a surface temperature of 280K out 
to 1.5AU from an early-type M dwarf star and 10 AU from a G type star. Neglecting the ef- 
fects of clouds and of gaseous absorbers besides H 2 , the flux at the surface would be sufficient 
for photosynthesis by cyanobacteria (in the G star case) or anoxygenic phototrophs (in the M 
star case). We argue that primordial atmospheres of one to several hundred bars of fi2-He are 
possible, and use a model of hydrogen escape to show that such atmospheres are likely to persist 
further than 1.5 AU from M stars, and 2 AU from G stars, assuming these planets have protecting 
magnetic fields. We predict that the microlensing planet OGLE-05-390L could have retained a 
H2-He atmosphere and be habitable at ^2.6 AU from its host M star. 

Subject headings: astrobiology — planetary systems 
1. Introduction 



type stars ( Koch et al. 201dl : Borucki et al. 2011 ) 



In the circumstellar habitable zone (HZ), an 
Earth- like planet has s urface temp eratures per- 
missive of liquid water ( Hard Il979f) . The HZ is 
typically calculated assuming an H2O-CO2 green- 
house atmosphere. At its inner edge, elevated wa- 
ter vapor leads to a runaway greenhouse state; 
at the outer edge, CO2 condenses onto the sur- 
face. The HZ of the Sun is presently 0.95-2 AU 
but is muc h more compact aro und less-luminous 



M dwarfs (jKasting et all 11993). Most reported 



but these cannot be confirmed by Doppler. In 
contrast, the microlensing technique, in which a 
planet around a foreground gravitational lensing 
star breaks the circular symmetry of the geom- 
etry and produces a transient peak in the light 
curve of an amplified background star, is sensi- 
tive to planets with projected separations near the 
Einstein radius of their host star, e.g. ~3.5 AU 
for a solar ma ss. Surveys have fo und planets as 
small as 3M® (jBennett et al.ll2008l) and pr ojected 



Earth-size exoplanets are well starward of the HZ 
because the Doppler and transit detection tech- 
nique s are biased toward s planets on close-in or- 
bits ([Gaidos et al.l 120071 ). The Kepler mission 



separations of 0.7-5 AU (jSumi et ai]|2010l) . Most 



can detect Earth-size planets in the HZ of solar- 



microlensing-dctcctcd planets orbit M dwarf stars 
as these dominate the stellar population. The Ein- 
stein rad ius of these stars is well outside the clas- 
sical HZ (jSumi et al.ll2010l:lKubas et al.ll2010l ) but 
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coincides with the "ice line" where N eptune-mass 



or icy super-Earth plan ets may form (| Gould et al 
l2010l : iMann et al.ll2010h . 



The effective temperatures of microlensing 
planets around M stars will be below the conden- 
sation temperature of all gases except for H2 and 
He. These mono- and di-atomic gases lack the 
bending and rotational modes that impart fea- 
tures to the infrared absorption spectra of more 
complex molecules. At high pressure, however, 
collisions cause H2 molecules to possess transient 
dipolcs moments and induced absorption. The 
resulting gree nhouse effe c t can considerably warm 
the surface. IStevenson (1999) conjectured that 
'free-floating" planets with thick H2-He atmo- 
spheres might maintain clement conditions with- 
out benefit of stellar radiation. Around a star, 
an fi2-He atmosphere might admit visible light to 
the surface: these planets could, in principle, host 
photosynthetic life. 

We demonstrate that H2 and He alone can 
maintain surface temperatures above the freez- 
ing point of water on Earth- to super-Earth-mass 
planets beyond the HZ. In Sf3]we describe calcula- 
tions of a radiative-convective atmosphere and its 
transmission of light; results are given in Ej3j In Sf4] 
we estimate the amount of H2 gas that a planet 
may accrete from its nascent disk, and the amount 
lost back to space. Finally, we discuss the effects of 
other gases and a candidate H2 greenhouse planet 
detected by microlensing 

2. Calculation of radiative fluxes and sur- 
face temperature 

Surface temperatures were computed using a 
one-dimensional radiative-convective model. For 
pure H2 and H2-He atmospheres, the infrared 
opacity d u e to c ollision-induced absorption is from 
Borvsow ( 20021 ). (Effects of additional green- 
house gases are discussed in ij5J) Thermal infrared 
fluxes for any given temperature profile were com- 
puted using the integral form of the solution to 
the non-scattering Schwarzschild equations, while 
the planetary albedo and absorption profile of in- 
coming stellar radiation were computed using a 
two-stream scattering code including the effects 
of continuum absorption and Rayleigh scatter- 
ing. The incoming stellar flux was modeled as 
blackbodies at 3500 and 6000 K for M and G 



stars, respectively. The calculations were indepen- 
dently performed in wavenumber bands of width 
10 cm -1 over 10-16490 cm -1 . The atmosphere 
was assumed not to absorb at higher wavenum- 
bers (A < 606 nm), although Rayleigh scatter- 
ing of incoming stellar radiation was still consid- 
ered. Continuum absorption is a smooth function 
of wavenumber, obviating the need for a statistical 
model within each band. Convection was modeled 
by adjustment to the adiabat appropriate for a 
dry H2-He composition (but see fj5]). The adiabat 
was computed using the ideal gas equation of state 
with constant specific heat. 

In selected cases up to 10 bar surface pres- 
sure p s , the full radiative-convective model was 
time-stepped to equilibrium. The atmosphere 
was found to consist of a deep convective tro- 
posphere capped by a nearly isothermal strato- 
sphere. Absorption of incoming stellar radia- 
tion slightly warms the stratosphere, but does 
not create an Earth-like temperature inversion 
because the pressure-dependence of absorption 
means that most absorption of stellar radiation 
and heating occurs at low altitude or at the sur- 
face. The tropopause occurs at a pressure level 
of <170 mbar, producing an optically thin strato- 
sphere that has little effect on infrared cooling to 
space. For p s = 1 bar, the formation of the strato- 
sphere cools the surface by only 3 K relative to a 
calculation in which the entire temperature profile 
is adiabatic. The slight effect becomes even less 
consequential as p s and the infrare d optical thick- 
ness o f the atmosphere increases ()Pierrehumbertl 
( 201dh . Chapter 4). Therefore, we used the com- 
putationally efficient "all troposphere" approxi- 
mation, in which the temperature profile is set to 
the adiabat corresponding to an assumed surface 
temperature T g . 

The radiation code calculates the correspond- 
ing infrared flux to space OLR(T g ,p s ) (Outgo- 
ing Longwave Radiation), and planetary albedo 
a(T g ,p s ). T g is obtained from the balance 

OLR(T g , Ps ) = ±(l-a(T g ,p s ))L (1) 

where L is the stellar constant evaluated at the 
planet's orbit (the analogue of Earth's Solar con- 
stant Lq). 
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3. Climate model results 

3.1. Minimum surface pressure to main- 
tain surface liquid water 

We present results for a 3M m rocky planet wit h 
surface gravity g = 17 m s -2 (jSeager et al . 2007). 
Gravity enters the calculation of OLR only in the 
combination p 2 s /g] the results can be applied ap- 
proximately to other values of g by scaling p s . 
Scaled results will diverge slightly from the cor- 
rect values because Rayleigh scattering depends on 
Ps/g, whence albedo scales differently from OLR. 
Figure [1] (left panel) illustrates the energy balance 
for a pure H2 atmosphere. OLR computed with 
the all-troposphere approximation is shown as a 
function of p s for different values of T g , together 
with curves of absorbed stellar radiation for G and 
M star cases with stellar constants L = 40 W m~ 2 
and L = 80 W m~ 2 . Equilibria are represented 
by the intersections between OLR curves and ab- 
sorbed stellar radiation curves. The absorbed stel- 
lar radiation decreases more rapidly with p s for G 
stars than for M stars because the greater short- 
wave radiation in the former case leads to more 
Rayleigh scattering, whereas the latter is domi- 
nated by absorption. For L — 40 W m~ 2 (some- 
what less than Jupiter's insolation) p s = 10 bar 
is sufficient to maintain a surface temperature of 
280 K about a G star; around an M star, only 7 
bars is required. 

Figure [1] (right panel) also shows the minimum 
p s required to maintain T g = 280 K as a function 
of the orbital distance. We use 280 K as a crite- 
rion because an ocean world whose global mean 
temperature is too n ear freezing is likely to en - 
ter a "snowball" state ( Pierrehumbert et al.ll2011 ) . 
The required pressure is higher around M stars be- 
cause more greenhouse effect is needed to compen- 
sate for the lower luminosity. For an early M-type 
star with luminosity 1.3% of the Sun, 100 bars of 
H2 maintains liquid surface water out to 2.4 AU; 
around a G star, these conditions persist to 15 AU. 
Cases requiring > 20 bars are speculative because 
the maintenance of a deep convective troposphere, 
when so little stellar radiation reaches the sur- 
face, is sensitive to small admixtures of strongly 
shortwave-absorbing constituents. 

Our results are not appreciably affected by the 
formation of a stratosphere. The stratosphere in 
these atmospheres is essentially transparent to in- 



coming stellar radiation and heated mainly by ab- 
sorption of upwelling infrared radiation. As ex- 
pected for a gas with frequency-dependent ab- 
sorptivity, the stratospheric temperature is some- 
what (few K) below the grey-gas skin tempera- 
ture ?jf f/2 1 / 4 where •■ oTt f; = OLR = |(1 - a)L 
(|Pierrehumberd (|2010h . Chapter 4). 



3.2. Photosynthetic Active Radiation at 
the Surface 

Starlight is a potential energy source for life, 
so we consider whether the stellar flux reaching 
the surface is sufficient to sustain photosynthe- 
sis. Most terrestrial photosynthetic organisms use 
light in the A = 400 - 700 nm range. The flux 
to sustain half maximum growt h rate for many 
cyanobacteria is ~1 W m~ " 2 (|Carr fc Whittonl 
1982c iTilzen 119871 ) , and an anoxygenic green sul- 
fur bacterium can use a flux of 3 x 10 -3 Wm" 2 



(jManske et al.ll2005l ). Planets on more distant or- 
bits experience less stellar irradiation and require 
a thicker atmosphere to maintain surface liquid 
water, which reduce transmission of light to the 
surface. For a pure H2 atmosphere, transmission 
in photosynthetically active radiation (PAR) is 
limited by Rayleigh scattering and differs slightly 
for M vs. G star spectra. The main difference is 
that G stars have higher luminosity, requiring less 
atmosphere to maintain liquid water at a given 
orbital distance, and a higher proportion of their 
output is PAR. Figure [2] shows surface PAR at the 
substellar point as a function of orbital distance 
for a minimum pure H2 atmosphere to maintain 
280 K. Hydrogen greenhouse planets as far as 
10 AU from a G star or 1 AU from an M star 
could sustain cyanobacteria-like life; distances of 
tens of AU or several AU, respectively, still permit 
organisms like anoxygenic phototrophs. 

4. Origin and Loss of H-He Atmospheres 

A body whose Bondi radius (at which grav- 
itational potential energy equals gas enthalpy) 
exceeds its physical radius will accumulate gas 
from any surrounding disk. The critical mass 
is less than a lunar mass and protoplanets may 
acquire primordial atmospheres of H g and He, 
as w ell as gases released by impacts (jStevenson 
19821 ). Earth-mass planets will not experience 
runaway gas accretion and transformation into 
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gas giants ( Mizuno Il980h . The atmospheres oJ' 
accreting Earth-mass planets will probably be op- 
tically thick an d possess an outer radiative zone 
(jRafikovl l2006h . [f the gas opacity is supplied 
by heavy elements and is pressure-independent, 
the atmosphere mass is inver sely propo r tional 
to the accretion luminosity ( Stevenson! 1982 : 
Ikoma k Gendall200rilRafikovll200fih . Inthe "high 
accretion rate" case (planetesimals are strongly 
damped by gas drag), the surface pressure will 
be 1.2a 1 ' 5 (M p /M( S ) bar, where a is the orbital 
semimajor axis in AU and M p is the planet 
mass. In the "intermediate accretion rate" case 
(gas drag is unimportant), the surface pr essure 
will be 16a 2 (M p /M e ) bar (jRafikovl [20061) . At 
3 AU, the predicted timescales for these two 
phases of accretion are ~ 10 6 and 10 7 yr, respec- 
tively, and bracket the lifetime of gaseo us disks 
(|Haisch Jr. et all 120011 lEvans et al.ll2009h . Thus, 
primordial atmospheres with surface pressures of 
one to several hundred bars are plausible. Hydro- 
gen is also release d during the reaction of wate r 
with metallic iron ( ElkinsTanton k Seagerlliooih . 



Giant impacts may erode t he proto-atmosphere , 
but n ot completely remove it ([Genda k Abell2003 , 



20041 ) . although the mechanical pr operties of an 
ocean may amplify such effects (jGenda k Abel 
2005). Instead, the low molecular weight fi of 
atomic H makes it susceptible to escape to space 
when stellar extreme ultraviolet (EUV) radiation 
heats the upper atmosphere, and by stellar winds. 
We examine EUV-driven escape of atomic H from 
a planet of radius r p using a thermal model of an 
H- or H2-dominated upper atmosphere, assuming 
the escape rate is not limited by H2 dissociation. 
The upper boundary is at the exobase where H 
escapes, and the lower boundary is where eddy 
diffusivity becomes important (in essence, the ho- 
mopause) . At the homopau se, the number de nsity 
is taken to be n = 10 19 m -3 ( Yamanakall995l ). and 
the temperature is set to the stratosphere temper- 
atures predicted by the radiative-convective model 
(see Sj3|). The exobase occurs at an altitude h e 
where the cross-section for collisions with molec- 
ular hydrogen is ~ 1, i.e. at a column density 
N = 2.3 x 10 19 m~ 2 . The exobase temperature T e 
and Jeans parameter A = GM p fi/((r p + h e )kT e ) 
are calculated by balancing the globally-averaged 
EUV heating qEUV with heat lost to conduction, 
radiative cooling, and escape. At A > 2.8 we use 



the Jeans escape rate (in surface pressure per unit 
time): 



(GM j 



(r p + h e ) 



/' 



2nkT f 



N e (l + \)e 



(2) 



where N e is the number density at the exobase. 
For Jeans escape, we assume H is the escape com- 
ponent and that H2 is efficiently dissociated be- 
low the exobase; this maximizes the escape rate. 
When A < 2.8 , marking the trans ition to hydrody- 
namic escape ( Volkov et al.ll20l"ol) . we assume that 
escape is energy-limited: 



qEUV- 



h t f 



(3) 



where h t is the altitude where EUV is absorbed 
(the thermosphere), and the right-hand factor ac- 
counts for thermal expansion. 

All EUV (A < 912A) is assumed to be absorbed 
in a layer at N = 2.3 x 10 19 m~ 2 . Lyman-a 
radiation resonantly scatters at higher altitudes 
(N ~ 1 x 10 17 m~ 2 ) but, in the absence of an ab- 
sorber such as ato mic oxygen, will escape b efore 
being thermalized (|Murrav-Clav et ah 20091 ). 

For given exobase and homopause tempera- 
tures, the temperature profile is calculated down- 
ward to its maximum at the thermosphere and 
then to the homopause, assuming a mean /i of 
3.1 (appropriate for H2+He in a solar propor- 
tion). The exobase temperature is adjusted un- 
til the correct value of n at the homopause is 
reached. Cooling from the escape of H at the 
exobase, the adiabatic upward motion of gas in 
the atmosphere to maintain stead y-state, and the 
collis i on-induced op acity of H 2 (|Borvsow et al 



19971: lBorvsowll2002l) are included. We use a ther- 
mal cond uctivity k » .0027 T 73 W m" 1 K~ 



based on 

qEuv — 



Allison k Smith! (|1971I ) 



We adopt 
AU from a 

4.5 Gyr-old G star (|Watson et al.l Il98lh . This 



1.5 x 10~ 4 W 



at 1 



is the current solar incident E UV flux of 1.67 x 
10~ 3 W m~ 2 for A < 900A (jrlinteregger et al. 



19811 : iDel Zanna et all 12010) tim es a 30% heat 



g efficiency (|Watson et al . 1981). EUV is scaled 



to different ages and M stars assuming propor- 



tionality to soft X -ray flux (|Penz et al 
Penz k Micelall2008h . 



2008; 



At sufficiently high EUV flux, the atmosphere 
expands to the scale of the planet, gravity at the 
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exobase decreases, no steady-state solution to the 
temperature profile exists, and A < 2.8. We pre- 
sume this to correspond to hydrodynamic escape, 
fix the temperature structure at its last state con- 
sistent with Jeans escape, and use the energy- 
limited escape rate (Equation [3|) . Additional H 
may esca pe as a result of inter action with the stel- 
lar wind ( Lammer et al.ll2007t > . 

We calculate total atmospheric lost by 4.5 Gyr 
for planets around G and M dwarf stars (Figure 
[3]). Significant atmospheric loss usually ceases by 
—2 Gyr. Planets around M dwarfs, with their 
lower X-ray luminosities, can retain their H2-He 
atmospheres close to the star. However, these 
stars, especially those of mid- to late-M spec- 
tral type, may yary greatly in their EUV flux 



(jReiners fc BasrJl2008t Browning et alJl201Clh 



5. Discussion 

Effects of other gases: The addition of He to 
the pure H2 atmospheres considered above would 
warm the surface by somewhat less than the ef- 
fect of adding the same number of H 2 molecules, 
because the CIA coefficients for H2-He colli- 
sion are generally we aker than those for H2-H2 



( Borvsow et al 



1988), and because He-He col- 
lisions do not contribute. Addition of He also 
makes the adiabat steeper and the surface warmer 
relative to the pure H2 case, but for a 10% He 
concentration the extra warming is no more than 
1 K. 

The warming provided by the H2 greenhouse 
permits other greenhouse gases to be stable 
against complete condensation, notably CO2, 
H 2 0, NH 3 and CH 4 . For the distant orbits of 
interest in this Letter, the effective radiating tem- 
perature occurs in the upper troposphere, and is 
-60 K for a planet 2 AU from an M star. The 
low saturation vapor pressure at this temperature 
confines the gases to lower altitudes, where they 
add little to the greenhouse effect. Even at 100 K 
the saturation vapor pressures of CO2, H 2 0, NH 3 
are all <0.02 Pa. CH4, with a vapor pressure of 
0.338 bar at 100 K, is most likely to exert a modest 
additional warming effect. 

The condensible gases also affect climate through 
latent heat release, which makes the adiabat less 
steep and reduces the surface temperature. Be- 
cause of the heat capacity of the thick H2 atmo- 



sphere, the effect on the adiabat is slight, amount- 
ing to, for example, a 3 K surface cooling for an 
H2O saturated atmosphere in 5 bars of H2. It is 
much less for thicker atmospheres. Instead, the 
main effect of condensible constituents would be 
to form tropospheric clouds, which would increase 
the albedo and hence increase the H2 pressure 
needed to maintain habitability. A quantitative 
treatment of clouds is outside the scope of this 
Letter. 

Cold cases: Two "super-Earth" -mass planets 
detected by microlensing are plotted in Figure [3Jd. 
MOA-2007-BLG-192L orbits -0.7 AU from a very 
late-type M dwarf. It has an effective tempera- 
ture of 40-50 K, below the condensation temper- 
ature of all gases except H2 and He, although N2 
or CO might be volatilized at the su bstellar point 



of a synchronously rotating planet (jKubas et al 
2010h . We predict that any primordial hydro- 



gen atmosphere has been lost by EUV-driven hy- 
drodynamic escape carrying He and other light 
volatiles with it. In contrast, more massive OGLE- 
05-390L, which orbits -2.6 AU from a mid Mi- 
type star but has a simila r effective tempera- 
ture ( Ehrenreich et al. 20061) . may retain a pri- 



mordial hydrogen atmosphere. This planet could 
potentially sustain liquid water at its surface, 
and may represent potentially ocean-bearing plan- 
ets to be revealed by future microlensing surveys 
(iBeaulieu et alJl2008t iBennett et al.ll201C 



All planets smaller than Neptune detected by 
Kepler or Doppler orbit within the zone where 
hydrogen escape is predicted to be efficient. 
Gliese 581d (> 6M e , 0.2 AU) is in this cate- 
gory; though it may lie within the classic HZ 



gory; tnougn it may lie witnm tne classic tlL 
(IWordsworth et alj l20ld; Ivon Paris etail l201Clt 
Hu fc Dindl201ll) . Larger, transiting planets such 



as GJ 436b (22 M®) have retained a low molecular- 
weight envelope despite their proximity to their 
parent stars, perhaps due to their high gravity 
and migration from further out in the primordial 
nebula. Should the Kepler mission be extended 
(—6 yr), planets at 1.5-2 AU, the inner boundary 
where a hydrogen atmosphere is retained (Figure 
[3|), could be confirmed with 3 transits. Measur- 
ing mass would be at the limit of current Doppler 
technology (radial velocity amplitude < 1 m s _1 ), 
but spectra obtained during transits might reveal 
any H2-rich atmosphere. 

The effects of life: These habitable worlds may 
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nurture the seeds of their own destruction, on 
account of the chemical disequilibrium between 
an H2-dominated atmosphere and (we presume) 
a comparatively oxidizing silicate mantle. If the 
planet's tectonics supports CO2 outgassing as is 
the case for Earth, then methanogens could de- 
plete the H 2 atmosphere by combining it with CO2 
to produce CH4. The O2 produced by cyanobac- 
teria would similarly consume H2 and convert the 
atmosphere to organic carbon and water. Either 
case would eventually cause a freeze-out of the at- 
mosphere, in the absence of some biotic or abiotic 
process that regenerates H2. 

This research was supported by NASA grant 
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between EG, RP, Scott Gaudi, and Sara Seager 
at the Aspen Institute for Physics. We acknowl- 
edge two anonymous reviewers of this manuscript. 
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Fig. 1. — Left panel: Determination of surface 
temperature from the top-of- atmosphere radiation 
budget. Solid curves show the OLR (infrared 
emission to space) as a function of surface pres- 
sure for the various surface temperatures indi- 
cated on the curves. Pairs of dashed lines give 
the absorbed solar radiation for stellar constant 40 
Wm -2 (short dashes) and 80 Wm~ 2 (long dashes). 
The upper curve in each pair is for an M star spec- 
trum while the lower is for a G star. Right panel: 
Surface pressure required to maintain 280K sur- 
face temperature, as a function of radius of a cir- 
cular orbit. Results are given for an M star (0.013 
times solar luminosity), and a G star (solar lu- 
minosity). All calculations were carried out for 
a pure H2 atmosphere on a planet with surface 
gravity 17 ms~ 2 Planetary albedos were computed 
assuming zero surface albedo. 
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Fig. 3. — EUV-driven atmospheric escape (bars) 
at an age of 4.5 Gry as a function of planet mass 
and distance from a G star (left) or M star (right). 
The heavy solid line is where mass loss equals 
the estimated mass of proto-atmosphere acquired 
by a planet during its early phase of oligarchic 
growth (high rate of planetesimal accretion) and 
the heavy dashed line is the same for the case 
of the late phase of oligarchic growth (interme- 
diate rate of planetesimal accretion). See text for 
details. Two microlensing-detected super-Earth- 
mass planets are plotted, with the boxes repre- 
senting the undercertainties in their parameters. 



Fig. 2. — Flux in the photosynthetically active 
band (400-700nm) reaching the surface at the sub- 
stellar point as a function of distance from a G or 
M host star. The flux is calculated assuming an 
H2 atmosphere with surface pressure sufficient to 
maintain 280 K surface temperature at each dis- 
tance. 
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